NPS ARCHIVE 
1969 

FRANSON, A. 

PREDICTION OF STATISTICAL 
SYSTEM PERFORMANCE FROM 
PARAMETER DISTRIBUTIONS 



by 



Alvin LaVerne Franson 







- 








'ATE school 

(3-5101 



United States 
Naval Postgraduate School 



PREDICTION OF STATISTICAL 
SYSTEM PERFORMANCE FROM 
PARAMETER DISTRIBUTIONS 



by 

Alvin LaVerne Franson 



June 1969 



Thti document kcu> been approved ^on pubtic tie- 
lexue and 6 ate; tt6 dti&Ubution t6 unlimited. 





SIS 




Prediction of Statistical System Performance 

from 

Parameter Distributions 



by 



Alvin LaVerne JFranson 
Lieutenant, United States Navy 
B, S. E. E. , University of Kansas, 1963 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN ELECTRICAL ENGINEERING 



from the 

NAVAL POSTGRADUATE SCHOOL 
June 1969 



c, 



C> 1 



N?S Av^H!V£ 

i <n 

•p^AKJSoK). A* 

ABSTRA CT 



Techniques for extrapolating statistical data on the parameters 
of individual components to the statistical performance index for an 
overall system are considered. Two cases are evaluated. (1) The 
deterministic case in which the system's performance index is known 
functionally in terms of the system parameters. (2) The non- 
dete rministic case in which only limited data on the performance index 
and its sensitivity with respect to system parameters is known. Com- 
puter programs are developed in both cases for combining given 
probability density distributions for the parameters into an overall 
probability density distribution for the system performance index. 
Theory and programs are developed and verified with a specific numer- 
ical example. 
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I. INTRODUCTION 



The performance of a system, from a statistical point of view, 
depends upon the nominal values of the parameters of the component 
parts included in the system and on the statistical distributions of 
those parameters. Alternately, the response of the component parts 
to a given signal may be known statistically and it is desired to pre- 
dict an overall statistical distribution for the system response. 
Statistical information on the system’s response may be interpreted 
in terms of probability of success or failure (operation or non- 
operation) for a single system or in terms of the percentage of 
operational systems that can be expected in an ensemble of systems. 
Currently three basic methods are being used to obtain performance 
data for systems; the worst-case method [Mark and Stember 1965], 
the moment method [Mark and Stember 1965], and the Monte Carlo 
method [Shreider 1966]. 

The worst- case method is a non- statistical approach intended to 
determine whether it is possible, with given parameter tolerance 
limits, for the systems performance characteristics to fall outside 
specifications. This method determines nothing of the probability 
distribution of the performance. Worst-case analysis, although it 
may be involved, does provide useful engineering insight and is used 
a great deal in practice. However, its outlook is basically limited. 
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Statistics are combined with system analysis techniques in the 
moment method to estimate the probability that performance will re- 
main within specified limits. The actual parameter probability 
distributions are usually replaced by Gaussian distributions having 
the same variances. Hence, the moment method has inherent limita- 
tions. Parameter probability distribution characteristics such as 
skewness and peakedness are usually not taken into account. 

In the Monte Carlo method a large number of alternate replicas 
of a system are simulated by mathematical models. Component 
values are randomly selected and the performance of each replica is 
determined for its particular set of components. The accuracy of the 
Monte Carlo method is dependent upon taking a sufficiently large 
number of replicas and can be both lengthy and costly. Monte Carlo 
methods are also used extensively in practice although they provide 
limited engineering information. 

The purpose of this thesis is to present a numerical technique 
for combining statistical data on several independent components into 
a single distribution function for the system performance. This is 
accomplished by considering two general cases. The deterministic 
case, where the functional relationships between the system parameters 
and the output are known functionally; and the non-dete rministic case, 
where these functional relationships are unknown but measured dis- 
crete data is available. In both cases the actual parameter prob- 
ability distributions are used in conjunction with numerical integration 
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techniques to obtain the system distribution function. One major 
advantage of the approach is that the known parameter distributions 
are quite arbitrary and may be histograms obtained from the results 
of experiments on individual components. 
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II. DETERMINISTIC CASE 



The object of the deterministic case is to investigate the dis- 
tribution of a system output when the output has a known mathematical 
relationship to several internal parameters. The basic assumptions 
are that all of the internal parameters are independent and that the 
probability distributions of the internal parameters, although arbi- 
trary, are known mathematically or given in terms of histograms. 

To commence the investigation, the general equations for various 
combinations of parameter distributions which may be expected in a 
typical system performance index must be obtained. These will 
basically consist of operations on the distributions of random variables 
which are interrelated by addition, subtraction, multiplication, and 
division. The basic mathematics for these operations are well known 
in the literature, at least in general form [Beckmann 1967]. 

Consider two independent random variables: X, assuming values 

x with a known probability density p^(x), and Y, assuming values y 
with a known probability density py(y)* Also consider the functions 
U = U(X, Y) and V = V(X, Y) and their inverse functions X = X(U, V) 
and Y = Y(U, V). Then, when X is near x and Y is near y, U and V 
must be near u = u(x, y) and v = v(x, y). Hence, 

P XY ( x > y) dxdy = p uv (u, v) dudv (1) 

Substituting the inverse functions and the required probability 
densities, we obtain 
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P uv (u, = p x f X ^ U ’ p Y f y ^ u ’ 



( 2 ) 




where Puy( u > v) is the joint probability density of U and V, p-^[x(u, v)] 
is the probability density of X, py [y(u, v)] is the probability density 
of Y, x(u, v) and y(u, v) are the inverse functions, and 
is the Jacobian, given by 



d< 


Jx,y) 


3< 


tu,v) 




ax 


ax 


au 


av 


11 


11 


au 


av 



(3) 



When X and Y are related by a single functional the foregoing can be 
modified by letting U = U(X, Y) and V = X. (4) 

Then, 




ax 

au 

11 

au 




Equation (2) then becomes 

p UX (u ’ x) = p X (x)p Y f y(u ’ x) ^|'U| 



(5) 



( 6 ) 



Since, 

as 



in general, py(u) 



■ J" p ux (u> 



x)dx, 



equation (6) may be written 



p U (u) = J Px (x)p Y t y(u ’ I'UI dx (7) 

- oO 

Using equation (7), the specific formulas of Table I have been 
compiled for the functional relationships and operations listed, which 
were chosen on the basis of what might be expected in engineering 
problems. 
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The technique outlined here consists of the successive application 
of the formulas of Table I as required, with the required integration 
performed numerically. Various types of numerical integration are 
possible; rectangular integration, trapezoidal integration, or Simpson's 
Rule. These are summarized as follows; 

Rectangular integration [Stanton 1961], 



r 

J : ' ( x ) s x : fjx. ) 

X, 



Trapezoidal integration [Stiefel 1963], 



( 8 ) 



f X " f ''j 

J f(x)dx Axj| [f(y. 1 )+f'(7 n )J +) f(>- i )l 

*> v i^-2 J 

Simpson's Rule [Stanton 1961], 



(9) 



J 



X 

f(x)dx 



-¥X>v 



where C 



.1 0 

1, 4, 2, --- 2, 4, 1 



( 10 ) 



Although the specific technique of numerical integration is of 
prime importance in so far as accuracy and computer programing 
and time are concerned, the exact procedure used is a detail which 
is best dictated by the particular problem to which these techniques 
are to be applied. 

This method may be best illustrated by considering the following 
example. Consider that: 



h = ^ 4 £ 



(in 
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where a and b are constants and X, Y, and Z are independent random 
variables with known (measured) probability densities. Consider that 
Px(x), Py(y)> and Py(z) are given as follows: 




(a) Consider first the combination U = XY, with = 0. 1. 

Using rectangular integration for simplicity, equation (3) of Table I 
becomes, for u = U-^ : 



When u 



Pu (u l) = Jp x ( x )py( u / X )^ 

-00 

Uj = 0. 09 



y^-Py (Xj,) 

i=l 



x. 



( 12 ) 



y>i 


U, /X; 


&(**) 

Xl 


ft ( a 'hi) 




0. 1 


0. 9 


4. 0 


1.8 


7. 2 


0. 2 


0. 45 


4. 0 


0. 9 


3. 6 


0. 3 


0. 3 


4. 0 


0. 6 


2. 4 


0. 4 


0. 22 


4. 0 


0. 43 


1. 72 


0. 5 


0. 18 


4. 0 


0. 37 


1. 48 


0. 6 


0, 15 


2. 6 


0. 30 


0. 78 


0. 7 


0. 13 


1. 7 


0. 26 


0. 44 


0. 8 


0. 11 


1. 0 


0. 21 


0. 21 


0. 9 


0. 1 


0. 44 


0. 20 


0. 088 



Pu ( Ul ) = (0. 1 )( 1 7. 918) = 1. 79 



sum = 17.918 
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when u r u^= 0. 2 



Xi 


Ug A; 


X; 






0. 1 


2. 0 


4. 0 


0. 0 


0. 0 


0. 2 


1. 0 


4. 0 


2. 0 


8. 0 


0. 3 


0. 7 


4. 0 


1.4 


5. 6 


0. 4 


0. 5 


4. 0 


1. 0 


4. 0 


0. 5 


0. 4 


4. 0 


0. 8 


3. 2 


0. 6 


0. 33 


2. 6 


0. 67 


1. 74 


0. 7 


0. 3 


1. 7 


0. 6 


1. 02 


0. 8 


0. 25 


1. 0 


0. 5 


0. 5 


0. 9 


0. 22 


0. 44 


0. 42 


0. 185 



Pu (u 2 ) = (0. 1) (24. 245) = 2. 4 



sum - 24. 245 



The pattern for evaluating p^j(u) is repeated for u = . 3, .4, 1. 0 

yielding the curve of Figure 1. 

... 1 l 

Here equation 

v - / j\. i 

(7) of Table I may be applied. 



(b) Next consider the function V = — = J— 

U XY 



v ; 



-p, (i v) . p '(u,l 



(13) 



jj 


v 


.'u > 
i i - * 


_ _ 


> ( - ’ 
V ’ • 


0. 1 


10. 0 


1.8 


0. 01 


0. 018 


0. 2 


5. 0 


2. 42 


0. 04 


0. 097 


0. 3 


3. 33 


2. 39 


0. 09 


0. 215 


0. 4 


2. 5 


2. 11 


0. 16 


0. 338 


0. 5 


2. 0 


1. 54 


0. 25 


0. 385 


0. 6 


1. 66 


1. 02 


0. 36 


0. 364 


0. 7 


1. 43 


0. 55 


0. 49 


0. 270 


0. 8 


1. 25 


0. 28 


0. 64 


0. 179 


0. 9 


1. 11 


0. 05 


0. 81 


0. 040 



Py(v) is plotted in Figure 2. 
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(c) The function W may now be calculated using equation (8) 



of Table I, where W 



a V + b 



a 

XY 



+ b 



( 14 ) 



and 





|a|M V ' 



If a = 2 and b = 1, the following result is obtained. 



(15) 



V 


W = aV + b 


Py( v ) 


P y( v) /2 


1 


3 


0. 0 


0. 0 


2 


5 


0. 385 


0. 1925 


3 


7 


0. 25 


0. 125 


4 


9 


0. 15 


0. 075 


5 


11 


0. 10 


0. 050 


6 


13 


0. 06 


0. 030 


7 


15 


0. 04 


0. 020 


8 


17 


0. 03 


0. 015 


9 


19 


0. 02 


0. 010 


10 


21 


0. 01 


0. 005 



The result is plotted in Figure 3. 

(d) Since, from equations (11) and (14), 

H = W + Z (16) 

it follows from equation ( 1 ) of Table I that pjj(h) is given by the con- 
volution of p^y(w) and p^(z). This operation may also be carried out 
numerically and for the given p^(z) may be expressed as 




Vt-I 



These integrals may also be carried out numerically yielding the 



results plotted in Figure 4. 
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To work more complicated problems than the previous example 
requires a great amount of mathematical calculations. Thus, when 
using numerical integration techniques the method lends itself nicely 
to the digital computer. Figures 5 through 9 are the computer flow 
graphs corresponding to some of the equations of Table I. The com- 
puter program for any mathematical function can then be derived 
from the appropriate combinations of these basic flow graphs. 

To illustrate a more complex problem consider the following 
example- 

given 

0 = 1285 f X * i) + W (18) 

where 

A = B 3a E 6a , B = B 3b B 6b 

and 



B 



3a 



B 3b 

1 + B 3b (2. 16 x 10-2) 



B 



B 



6a 



6b 



1 + B 6b (0. 63 x 10-2J 



is a normal distribution; mean = 130, standard deviation = 10 
3 b 

B^ is a normal distribution; mean = 65, standard deviation = 12 
W is a uniform distribution from w = - 4/7 to vv = + 4/7 
Theta is system output parameter and the variables B^^, B^, and 
W are independent. 
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Since the variables are independent, equation (18) can first be 



simplified by ordinary mathematical manipulation. 

Let X = B 3b , Y = B 6b , k } = 2. 16 x 10" 2 , k 2 = 0. 63 x 10“ 2 

then B = XY 

A = XY 

(l + k x X) (1 + k 2 Y) 



1 1 _ 
A " B 



k k + + _^2 

Y X 



and 9 = 0. 175 + 11 : A + JLA + w 

Y X 



(19) 



Now, utilizing the basic methods and flow graphs, the problem 
can be programed on the digital computer and the probability density 
of 9 readily obtained. The program and solution are contained in 
Appendix I. 
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III. NON- DETERMINISTIC CASE 



In many practical problems the functional relationship between 
variables is unknown and it cannot be expressed simply by formula. 

It is generally available only as numerical data from a computer run 
on a mathematical model of the system under study. Two cases are 
considered here. First, the case where the system response of 
interest, U(X), and its sensitivity coefficients, U’(X), are available 
at discrete values of a single variable, X. 



Second, the case where the system response of interest, U(X, Y), 
and its sensitivity coefficients, U^X, Y), are available at discrete 
values of two independent variables, X and Y. 



U(X) 




( 20 ) 




u 




( 21 ) 



U 11 u 12 U ln 



U(X. Y) 




( 22 ) 




u 



mn 



where 



u ij = u ( x i’yj) 



and 



20 



U'(X. Y) = 



u 



where 



11 


u 12 ‘ ' 


- u In 


21 


u< 22 “ - 


- u' 

2n 




- 


- 




- 


- 


ml 


U 'm2 


- u l 

mr 




au(xj 


■ y j ) 



( 23 ) 



dx i 

In both cases the problem reduces to determining p^(u) from the 
known distributions p-^(x), p^(y) and the data of equations (20) and (21), 
or (22) and (23). 

In the first case, consider that the probability density, p^(x), 
and the cumulative distribution. 



/ *j 

P x (x)dx 



(24) 



are known. Then at each data point of u(x). 



U 



[ u(x i'] 



= P x (Xj) 



(25) 



That is, the probability that u — u(Xj) * s e< l ua ^ the probability that 
x ^ Xj . Equation (2 5) gives discrete points on the cumulative 
distribution, P^j(u), for each value of the parameter x for which 
u(x) has been measured. 

It is also well known [Beckman 1967] that 

Py \ - / 



Pu(u) = 






' AU 


1 ,A~i 1 




1 d x 1 



(26) 
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Thus the slope of the cumulative distribution curve. 



\ * ( : 
Jl ^ I 

ov . 



at each data point is equal to the known probability density, p-^(Xj), 



CV 

f' V 



, at the 



divided by the measured sensitivity function, 
data point. The conclusion is that the cumulative distribution, Py(u), 
and its slope may be calculated readily from the measured response 
and sensitivity data. This curve is therefore evaluated as a piece- 
wise linear approximation. 

When U is a function of two variables, X and Y, the foregoing 
cumulative distribution curve can be evaluated for each value of the 
second variable, Y. Thus a family of such curves or data points 
denoted by Py(u, yj) is available. For a specific value of u = u. 
each of these curves gives the probability that u ^ u^ when y = y y 
The total probability that u ^ m depends upon Py(y)Ay and is 



given by the equation: 



m 



Pu(ui) = P u (u i ,y H )p v (vJ 

. 1=1 



(27) 




constant and x as a single variable as discussed previously. 



This result can be proven by considering the following (see 
equation (7) ): 



Pyfu) = 



r 

U, 



'.v)p x [-/.(u,y)] | 



( 28 ) 



J 

- oO 

where u = u(x, y). ^ is the sensitivity function with respect to 

o 

the parameter x and in equation (28) must be expressed in terms of 
u and y. It is convenient to denote 
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dU / X 



( 29 ) 



and then equation (28) becomes 

CO 



Pu(u) - Jp Y (i')p x [x(l..y)] s ~^y T 



(30) 



We also know 

Py (U) = rp,r(u.)dU 

-00 

Substituting equation (30) into equation (31) and inte r- changing the 
order of integration yields 



(31) 



p u< 



u) ' f^(y)ayfpx[ x < u -y)]efSorT 

J J A 

— 00 — 00 



(32) 



Now 



r 


1 


F x|_ 


£(u,y)J 



Pg(u,y)du -- P Ify (u,y) (33) 



- 00 



- OO 



Since, if Y is considered a constant and X as the random variable, 
u = u(x, y), 



pu( u ’ y) = 



Px [ x ( x > y)] _ Px N*. 



(34) 



Ah 

dx 



u x (u, y) 

Substituting equation (33) into equation (32) yields: 



p u( u > = 



r 

I Py(y)P U (u, y)dy 

- 0O 



(35) 



In discrete form, equation (35) becomes equation (27) when 
rectangular integration is used* However, more accurate integration 
techniques may also be applied to equation (3 5), 
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To demonstrate the application of this non-dete rministic method 
consider part of the previous example of equation (18). It is now 
assumed that the functional relationship between the output and the 
parameters is unknown. However, we do know the probability den- 
sities of the parameters x and y and have certain measured data which 
consists of measured values of x, y, and 0 and the corresponding 
sensitivities for each measured 0. The data is presented in 
Table II. Following the procedure outlined, the problem can be 
programed on the digital computer and the probability density of 0 
readily obtained. The program and solution are contained in 
Appendix II. The results can be compared with those obtained using 
the deterministic formulas. Noting, of course, that the solution 
to the deterministic problem } contained in the Appendix, includes 
the addition of the uniform distribution W(w) and the solution to 
the non-deterministic problem does not. 
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IV, CONCLUSIONS 



The numerical technique presented here has distinct advantages 
over the worst-case, the moment, and the Monte Carlo methods. By- 
utilizing the actual parameter probability distributions the results 
obtained are more accurate than those of the moment method and pre- 
sent a great deal more information than the worst-case method. The 
accuracy obtained is equivalent to that of a large number of Monte 
Carlo replicas. Thus the saving in time over the Monte Carlo method 
is obvious. In general, the numerical method presented here can 
obtain a more accurate solution in less time than any of the other 
current methods. However, certain limitations are inherent in the 
basic assumptions used to formulate the method. Since independent 
parameters are required, a great number of practical problems can 
not be handled by this method as it is presented here. By extending 
the mathematical derivations to include the correlation of parameters, 
the method could be used for systems in which the parameters are 
not independent. 
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Table I 



(1) U 

(2) U 

(3) U 

(4) U 

(5) U 

(6) U 

(7) U 

(8) U 



Functional Formulas 



= X + Y 


c °* 

PyCu) = j p x (x)p y (u-x)dx 

• so 

= Jp x (u-y)p Y (y)dy 

• «o 


= X - Y 


r 00 

p y (u) = J p x (x)p y (x-u)dx 

- <o 

= J'p x (u+y)p y (y)dy 

. so 


= XY 


Pu (u > - J 

- OC 


= X/Y 


1 r°° 

p y (u) = — ^ j p x (x)p y (x/u) Ixldx 

= J P x (uy)pY(y) lyldy 

. oo 


= X 2 


Py( /U) + P y (-/U) 

P u (u) - u - 0 

U 2/u 


= IXl 


Py(u) = P x (x) + P x (-x) X-Q 


= a/X 


Py(u) = —^2 P x (a/u) 


= aX + b 


Pu (u) ■ -I- Px(^) 



(a and b are constants) 
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Table II 



Sample Data 



>5J 


41 


53 


65 


77 


89 


110 


0.9266 


0.7731 


0.6761 


0.6096 


0.5616 


120 


0.9204 


0.7669 


0.6699 


0.6034 


0.5554 


130 


0.9152 


0.7617 


0.6647 


0.5982 


0.5502 


140 


0.9108 


0.7573 


0.6603 


0.5938 


0.5458 


150 


0.9069 


0.7534 


0.6564 


0.5899 


0.5419 


a 


41 


< 

53 


> 

65 


77 


89 


*0 


0.0165 


0.00988 


0.00658 


0.00469 


0.00352 



X 


110 


120 


130 


140 


150 


<>x 


0.000670 


0.000561 


0.000478 


0.000413 


0.000359 
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3 



PjjU) 




Figure 1 




Figure 2 
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Figure 4 
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U = a/X 





Calculation of u^ and where U = a/X and a, x 

gmd p(x i ) are known. 



Figure 5 
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To perform rectangular integration of the distribution p(x^) 
vs. the distribution must be approximated by rectangular 
areas of equal widths ax. This program calculates the x^'s 
such that ~ x j + and interpolates for the correspond- 

ing values p(x^) from the given distribution p(x^) vs. x^ 



Figure 6 
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U = X + Y 




32 




This program performs the convolution U = X + Y by using 
graphical convolution techniques. The distributions p(x^) 
vs. x i and pCy^) vs. y^ are represented by sets of rec- 
tangular areas of width ax. The distribution of x^ is 
reversed and passed across the y , distribution in ax in- 
crements. The distribution vs. u^ is then obtained 

by the normal graphical convolution mathematical calcula- 
tions [Bracewell 1965 > Healy 1969J • 



Figure 7 
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U = XY 




34 




Figure 8 
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U = X/I 
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Figure 9 
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APPENDIX I 



The solution to the deterministic problem has been obtained by 
using three separate programs. Linear interpolation and rectangular 
approximations have been used to perform all convolution and rec- 
tangular integration. 

The problem and given information is as follows: 

Obtain the cumulative distribution of 0 where 



a _ p j_ A C . 

C = 0.175 
A = 8. 1 
B = 27. 8 
X(x) is N( 130; 10) 

Y(y) is N(65;12) 



tf(w) 



W(w) is uniform 0. 875 from - 4/7 to + 4/7 



Program one inverts the distribution x vs. p(x) and multiplies 
it by the constant A. Program two does the same for B/Y(y). Since 
both X(x) and Y(y) are normal distributions, the input data for these 
programs was taken from standard normal tables. 25 data points 
were used in each, the median and each 1/4 standard deviation out to 

t i 

3 standard deviations. The resulting distributions, X (x) and Y (y), 
are then used as input data for program three. Program three per- 

i t t 

forms the convolution U(u) = X (x) + Y (y), the addition U (u) = U(u)+ C, 
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and the convolution 0 = u’(u) + W(w) to obtain the cumulative dis- 

tribution of 9. 

To obtain highly accurate results a of 0. 002 was used for 

convolution and rectangular integration. The over-all accuracy for 
this type of solution will depend upon how closely the rectangular 
approximations resemble the continuous distributions which they 
represent. 
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A P P P N 0 1 X I 



PROGRAM ONF 



C FVAIUATION OF THE FROBAGILITY DENSITY FUNCTION, P(X'), 

c WHOKF X ' = A / X AND A, X, AMO P ( X J ARE KNOWN. ? ( X ) HAS 

C N( P b't L DISTRIBUTION. N = NUMBER OF POINTS X(I), m IS THE 

C. STANDARD DEVIATION OF THE DISTRIBUTION P(X). 

DIMENSION X ( 50) ,PX(50) ,U< 50 ) ,PU( 50) 

R E A T ( 5 , 1 C 0 ) N 
R S AD( 5, 1C1 )M. 

Pf ADI 5, 102) A 
100 FCRVATI 12) 

10 I FCR MAT (I 2) 

102 FPPMATIF 10. 0) 

WRITi(f ,210) N 

210 F 0 R M A T ( 1 H 1 , T 2 0 , ' N = *,I2//) 

WP I TF < o , 2 1 1 ) M 

211 F 0 RMA T ( T2 C , ’ M = ',12//) 

WRITt (6,212) A 

212 FrpMAT(T20, 'A = ',F10.5//) 

WRITE (6,213) 

213 FCPMAT(T20, * X ( I ) ' ,T50, 'PX( I ) ' / / ) 

DO 300 I = 1 , N 

° r / D ( 5 , I C 3 ) X { I > , p X C I ) 

103 Ff RMA T ( 2 F 1 1 .0) 

WP ITF (6 ,2\n X ( I ) , PX ( I ) 

2 1 ^ F0FMAT(2F 3C. B) 

3 CO CONTINUE 

DO 301 1 = 1 , M 
K=N-I +1 

F U ( K ) = ( X ( I )**2/A)*(PX( I ) / M ) 

U ( K ) = A / X ( I ) 

301 CONTINUE 

WP ITE <6 ,2C0) 

2 0"' FORMAT! //,T20 ,' X' ,T50, ' P( X )',//) 

DC 302 1=1, N 
201 FORMAT! 2F 3°. ft ) 

STOF 
[.! D 

N = 2 f 
M = 10 
A = F. 1 



X ' 



PX* 



C. 50624">96E-01 
'". 514285685-01 
0. 52258059E-01 
C. C 3114749E-01 
C.53999994E-01 
C. 54Q15249E-01 
0. 5 5 B62 C625-0 1 
C • 56 842 1 00F-01 
C.57R57137F-01 
C. 58509085E-01 
C. 59999995F-01 
C.61 132070E-01 
C. 62307686F-C1 
0. 63529372E-01 
C. 6A799964F-01 
0.581224135-01 
C. 674999955-01 
C . 68 C 3611CE-01 
0. 7043^-749R-01 
(,. 71Q9R967E-01 
0. 736363532-01 
0. 75348794E-G1 
0. 77142835E-01 
C. 79024374E-01 
C. 809999700-01 



0. 1 3 °04 1 640 01 
0.2786»73b0 01 
0.519058610 01 
0.9101 51 96E 01 
0.149999956 0? 
C.23179779E 02 
0.33614029F 02 
0. 457747940 02 
0. 5355 RQ 140 02 
0.70279877E 02 
0.79222473E 02 
0.833148190 02 
0.832272640 02 
0.776085360 02 

C. 679205170 02 
0.55782471E 02 

D. 4J0222C2F 02 
0. 31 123 7 03E 02 
0. 21 143661= 02 
C.134P4363E 02 
0. 80666637E 01 
0 . 45226297E 01 
0. 23319437F 01 
0. 11803312 F: 01 
0 . 5 4 3 2 0 9 6 7 F 00 
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noon 



APPENDIX I 



PROGRAM TWO 



EVALUATION IF THE PRCBA8I L I TY DENSITY FUNCTICN, P(Y')* 
WHERE Y'=B/Y AND 3, Y, AND PIY) ARE KNOWN. ? ( Y ) HAS 
NORMAL DISTRIBUTION. N=NUMB£R OF POINTS Y(I), M IS THE 
STANDARD DEVIATION OF THE DISTRIBUTION P(Y). 

DIMENSION X(50)»PX(50) ,U(50) ,PU(50) 

REALM 5, IOC )N 
Rt AD! 5, 101 >M 
Rt AD ( 5 , 1 C2 ) A 
100 FORMAT! I 2) 

L 0 1 FORMAT! I 2) 

102 FDRMAT(FIC.O) 

WRITE! 6, 210) N 

210 FORMAT! IHI ,T20 , *N = ',12//) 

WRITE (6*211 ) M 

211 FORMAT ( T 20 , ' M = • ,12//) 

WRITE (6,212) A 

212 FORMAT ( T20 , ' A = «,F10.5//) 

WRITE (0,213) 

213 FORMAT! T20, • Y ( I ) • ,T50, 'PY( I ) •// ) 

DO 300 1 =1 , N 

READ!?, 103) X(I),PX(I) 

103 FORMAT! 2F1C.0) 

WP I TE ( 6 , 2 l z ) X! I > ,PX( I ) 

2 1 A F0PMAT(2E3C.8) 

300 continue 

DC 3 C 1 I = 1 , N 
K = N— I + 1 

PU(K )=( X! I ) **2 / A ) *( PX( I ) /M) 

U(K)=A/X( I > 

301 CONTINUE 
WRITE (6,200) 

200 FORMAT!//, T20, 'Y* ,T50, 'P(Y )',//> 

DC 302 1=1 , N 

201 FOR MAT! 2E3 C. 8 ) 

STC1F 

END 
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27. B 








Y • 




PY • 




G.2752474QE 


00 


0. 1345455 IE 


00 


0. 2 836732 9£ 


00 


C. 2619796<+E 


00 


0.29263145E 


00 


0.47343379F 


00 


C. 30217373E 


00 


0.60A29284E 


00 


0. 31235939E 


00 


0. 12821751E 


01 


0.32325566= 


00 


0. 19132919E 


01 


0.33493960F 


00 


0.26742373F 


01 


0. 34749985E 


00 


0.35031176E 


01 


0. 361C3874E 


CO 


0.43010139c 


01 


C.37567550E 


CO 


0.4^4251736 


01 


0.391 5491 1 E 


00 


0. 5320552BE 


01 


G. AQ882331E 


00 


0. 53600159F 


01 


0.42769212E 


00 


0.50520163E 


01 


0.44838685E 


00 


0.4455861 IF 


01 


0.47118622E 


00 


0.35740417E 


01 


C. 496428316 


00 


0. 2330484*= 


01 


C. 52452803c 


00 


0.2 0377035E 


01 


0. 55599976E 


00 


0. 13694053E 


01 


C. 59148908E 


00 


0. 8575 1051E 


00 


C. 63181788F 


00 


0. 50082946E 


00 


C. 67804843E 


00 


0. 27210426E 


00 


0. 73157859E 


00 


0. 13721460E 


00 


C. 79428536E 


00 


0. 64261079F' 


-01 


0. 86874962E 


00 


0.27932845E- 


-01 


0. 95862025E 


00 


0. 11092324F- 


-01 
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APPENDIX I 



PROGRAM THREE 



CALCULA TION OF U=X + Y WHERE X, P(X), Y, P(Y), ARE GIVEN 
AND THE GIVEN VALUES CF X ( I ) AND Y ( I ) ARE INCREMENTED 
BY DX. 

DIMENSION X( 500) , PX ( 5 CC ) , Y ( 500 ) , PY < 500 ) , U ( 530 ) , PU(500 > 
P W( 600 ) ,Z(1000) ,PZ(1000) » c ( 500 ) ,W(600) 

R f AD ( 5 * 1 0 1 ) OX 
131 Ff RMATIF10. 0) 

WR I TE ( 6 » 2 1 C ) OX 

210 FORMAT ( 1H1 , T2 0 , 'DX =',F10.3//) 

R E AD ( 5 » 1 C6 ) C 

106 FGRMAT(FIC.O) 

WRITE (6,213) C 

213 FCPMAT(T20, 'C =',F10.8//) 

P. ir A 0 ( 5 , 1 0 5 ) N 
105 FORMAT! I 5 ) 

WFITF <6,21 1 ) N 

211 FCRMAT(T20, 'N = • ,15//) 

R £ AD ( 5 , 1 04 ) m 

104 FORMAT! I 5 ) 

WFITF (6,212) M 

212 FORMAT! //,T20, *M = • ,15//) 

WRITE ( c , 2C2 ) 

2 02 FOPMAT( //,T20,'X!I)',T50,'PX(I)'//) 

CALL ONF ( X , PX) 

WF I TF ( f , 203 ) 

203 FCPMAT(//,T2 0,'Y(I)' ,T50,'PY(I)'//) 

CALL ONE ( Y , P Y ) 

K = N+ M- I 

U( 1 )=X( 1 ) + Y ( 1 ) 

DO 302 1 =2 , K 
302 U ( I ) =U! I - 1 ) +DX 
DO 3C3 J = 1 , K 
PU( J)=0.0 
DO 3CF 1=1 , I 

I F( ( J-I+l ) .GT.M) C,C TO 30 A 
IF ( I .GT.M) GO TC 303 
PU!J)=FU(J)+PX(I)*PY(J-I+1)*DX 

304 CONTI NUF 
30 3 CONTIMJF 

C CALCULATION OF L(J)=U(J)+C WH£PE 'C* IS A CONSTANT 

DO 3C6 I = 1 , K 
U( I )=U( I ) +C 
30 6 CONTINUE 

WRITE (6,200 

2 00 FORMAT! U'l ,T20, 'U(I) * , T50 , ' PU ( I )'//) 
nr 305 1 = 1 , k 

305 WR I TE ( 6 , ?C 1 ) U ( I ) ,PU( I ) 

201 FC pm/t( 2E 3C. 8) 

C CALCULATION AND PLOT CF CUMULATIVE D ISTR I BUT I CN 

F ( 1 )=PU< 1 ) *DX 

DO 320 I =2 , K 

F ( I ) = F ( 1-1 ) + Pll( I ) *QX 

320 CONTI N1JF 
WPITF! 0225 ) 

225 FORMAT! 1H1 ,T20, 'U(J) ' ,T40, 'CUMULATIVE DISTRIBUTION') 

DO 321 I =1 , K 

WRITE! 6,201 ) U ( I ) , F ( I ) 

321 CE NT I NUF 

PERFORM THE CONVOLUTION OF 'U' WITH THE CCNSTANT PD c 
•W' (Z=U+W). 

CALL TWC! W , P W , L ) 

NM=K+L-1 
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I ( 1)=U( 1 ) + W( 1 ) 

or 310 1=2 , nm 

310 Z < I )=Z( 1-1 ) +DX 
or 311 J = 1 , NM 
P7 ( J ) =0. 0 

or 312 I =1 , J 

I F ( ( J-T+1) .GT.L) GO TO 312 
I F ( I.OT.K ) GO TO 311 
P7(J)=P7(J)+PU(I)*PW(J-I+l)*OX 

312 CONTINUE 

311 CONTINUE 
WRITE (6,230) 

>30 FCPMAT( 1U1 ,T20 , *Z( I ) • ,T50, • PZ< I ) •// > 
nr 313 i=i , nm 

313 WPITF (6,231 ) Z( I ) ,PZ( I ) 

231 FCRMAT( 2E3C. 8) 

C C A LCUL " T I ON OF THE CUMULATIVE DISTRIBUTION OF 'Z*. 

P 7 ( 1 >=PZ( 1 ) *DX 

DO 33C 1=2, NM 

P 7 ( I )=F7( I-1)+PZ( I )*DX 

330 CONTINUE 
WPITF(6,220) 

220 FOPMAT( 1H1 ,T20, *U( 1 )• ,T40 , 'CUMULATIVE DISTRIBUTION') 
DO 331 1=1 , NM 

221 FORMAT ( 2E 3 C. 8 ) 

331 CONTTNUF 
c TOP 
END 



SUBROUTINE ONE ( X E , PXF ) 

C FF ON GIVEN VALUES CF X ( I ) AND PX(I), (THESE ARE U(I) 

AND PU(I) PREVIOUSLY CALCULATED), AND A CHOSEN OX, WE 
OBTAIN THE VALUES OF X(J>, WITH INT5VALS CX, AND 
CALCULATE THE CORRESPONDING VALUES OF °X(J) BY LINEAR 
INTERPOLATION. 

DIMENSION X( 500) ,PX(5C0) ,XP(500 ) ,PXE( 500 ) 

F LAD (5,100 N 
IOC eORMAT(I2) 

P TAD (5, 101) DX 

101 FORMA T ( F 1 0. 0 ) 

DO 3C0 1=1, N 

300 PF AD( 5, 1C? ) X( I ) , PX ( I ) 

102 FORMA T( 2F1 0. 0 ) 

1 = 1 

J=1 

Xf ( J)=X( T ) + 0 X 7 2 

400 I F ( XE : ( J)-X(I ) ) 500,501 ,5C2 

502 1=1+1 

I F ( I.LT. ( N+l) ) GO TO 400 
GO TO 203 

301 P XE ( J ) =PX( I ) 

GO TO 200 

5 CO PXE( J) = PX( 1-1 ) + (PX( I ) — P X ( I — 1 ) )/(X( I » — X ( I - 1 ) ) * ( X E ( J ) - 
X f ! - 1 ) ) 

2CC WR I TE ( 6 , 2C 1 ) XE(J), PXE(J) 

201 FORMAT ( 2F 30. 3 ) 

XF(J+1)=XE(J)+DX 
J = J+1 
GO TO ^00 
203 RETURN 
END 



SUBROUTINE TWO(W,PW,L) 

C THIS PROGRAM TAKES THE GIVEN CONSTANT PROBABILITY 

C FUNCTION 'W, WITH AMPLITUDE 0.R75 FROM -4/7 TO +-4/7, 
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AND OBTAINS THE VALUES CF W(I) WITH I NCRE M ENT CX=C.C02 



DIMENSION W(600), PWI600) 

W( 1 ) =-0.571 
FH 1 1=0. 875 
PX=C.002 
L = 57 1 

WP I T£ ( 6 , 2 50 ) 

2 5C F OP N’A T ( 1 H 1 , T20, • W (I) • ,T50 , • PW ( I ) • // ) 
nr 350 1=2, L 
W( T )=W( 1-1 l+DX 

rwu )=pw(i-n 

353 CONTINUE 

DP 351 1=1 , L 

WF ITE ( 6,251 > W(I) ,PW(I) 

51 FCPVAT(2E30. 8) 

51 CONTINUE 
K E TUP N 
f NO 



N = 15 

M = 545 

C = 6.175 
DX = 0.002 
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THETA 



0.631 28109E-Q1 
C. 65128075E-01 
C. 641280415-01 
0.621280456-01 
C • & Q 1 280445-0 1 
C . 5912RG44E-01 
0. 56 1289445-01 
0.541280435-01 
0 . ^2 1 2 8043E-01 
C. 501280435-01 
0.481280425-01 
C .^61280426-01 
0 . 44 12 8 042 E -0 1 
0.421280416-01 
0. 401280416-01 
C. 33123 04 IE-01 
C .361280406-01 
f . 34128040^-01 
C • 32 1 2 8040E-0 1 
C.301 28039F-01 
0. 29128039E-01 
0. 26128039E-01 
C.24128038E-01 
0. 22 1 28038E-01 
0. 20128038 E -01 
C. 181280375-01 
C. 16128937E-01 
0.141280376-01 
0.121280366-01 
C. 10128036E-C1 
C. 31230358F-02 
C • 61 280355E-C2 
■C . 41 2 803 52 6-C2 
0.212803535-C2 
C. 129035436-03 
0. 18719644E-02 
C.38719643E-02 
0.53719628E-02 
C. 7 3719594E-02 
C. 937195605-02 
C. 11871 953E-C1 
C. 13871 9496-01 
C. 15871946F-01 
0.178719426-01 
0. 19871 939F-01 
0.218719366-01 
C. 238719326-01 
C.25871929E-01 
C. 2 7 67 1 925F-01 
0.29871922E-01 
0.31 871 919E-01 
0.33971915E-C1 
C. 358719 1 2E-0 1 
0. 37871 908F-C1 
0. 398719056-01 
C. 418719026-01 
0.43871 3985-01 
0.45871 8955-01 
0.478718915-01 
0.498713885-01 
0.518718855-01 
0.53871 381E-01 
0. 55871 878E-01 
C. 578718 74F-0 1 
0. 598718715-01 
0.618718685-01 
0.638713615-01 
0.658713355-01 
0.678713095-01 



CUMULATIVE 
C I STR I PUT ION 

0.35 1622805-08 
0. 243545595-07 
0.93342241E-07 
0.2936691 3E-06 
0. 71189^305-06 
0. 14791O67E-05 
0.2734=2415-05 
0 • 462032685-05 
0.727990216-05 
0. 10S5R°82E-04 
0. 155035295-04 
0.213980275-04 
0. 296690C5E-04 
0.375373065-04 
0.4S194212E-04 
0.60356095E-04 
0. 757554956-04 
0. 931420655-04 
C. 113283375-03 
0. 13646450E-03 
0. 162986975-03 
0. 19 3167 89E-03 
0.2273402 IE-0 3 
0.265853015-03 
0. 33907127E-03 
0 . 35737362E-03 
0.411151216-03 
0.470806615-03 
0. 53675333 E- 03 
0.609415585-03 
0.63922667F-03 
0. 776627805-03 
0. 972066015-03 
0.975992535-03 
0. 109886165-02 
0. 121 11294E-02 
0. 13432535E-02 
0. 148569236-02 
0. 163S9030E-02 
0. 190333906-02 
0. 19794477E-02 
0.216766875-02 
0.236843415-02 
0.258216845-02 
0. 280928915-02 
0. 30502048E-02 
0 • 33053122E-02 
0. 35749928E-02 
0. 38596098E-02 
0.41595064E-02 
0.44750199E-02 
0.49064664E-02 
0.51 54 15 52E-02 
0.55193806E-02 
0. 539941445-02 
0.62975064E-02 
0.671 288 3 7E-02 
0.7145751 3E-02 
0.75962991E-02 
0.80647096E-02 
0.855114315-02 
0. 90557523E-02 
0.957866395-02 
0.101199865-01 
0. 10679796E-01 
0. 112581515-01 
0.118550925-01 
0.124706515-01 
0. 1 3 104849E-0 1 
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0. 59P71783E-01 
0. 71871753E-01 
0 • 73 8 71 732 E-0 1 
C.758717C6E-C1 
C. 778716805-01 
C. 79371 6556-01 
0.31 8716295-01 
C. 33871 603E-01 
C. 85671 577E-01 
C. 87871 552F-C1 
C. 69971 5265-01 
r.Ql 871500E-C1 
6. Q3871474E-01 
C. 958714495-01 
0.97871423F-01 
r.99«713976-0i 
C.10187137E 0C 
C. 103871356 00 
C.10537132E 00 
C. 10787129c 00 
C.10987127E 00 
C • 111871246 00 
0. 1 1387122E 00 
C.11587119E CC 
C.11787117E 00 
0.119971146 00 
0.12 1871 HE 0 3 
C.123971C9E CO 
C.12537106E CO 
C.127871C4E 00 
C. 12 G 671C1E 00 
0. 13187C99E 00 
0. 13337096E 00 
C.13587G93E 00 
0. 13787091E 00 
C. 13967088E CO 
C . 141 97036E 00 
C. 14387083E CO 
C.14587080E 00 
C. 14 7 97 078 E OC 
C.14387075E 00 
C.15187C73F 00 
0. 15387C7CE CO 
0. 15587068E 00 
C.15737065E OC 
C.15987062E 00 
C.1618706CE CO 
0. 1639705 7E 00 
0.165870555 CO 
C. 167970526 00 
C. 16987050E 00 
C.17187047F 00 
C . 17 387 044 E 00 
C. 17587042 F. 00 
0. 17787C39E 00 
0.17°87037E 00 
C. 18187034E 00 
C. 133870325 00 
0. 13587029E 00 
0. 187870266 00 
0. 18987024E 00 
0. 19187021E OC 
0. 1936701 9E 00 
C. 1959701 6E 00 
0.19787014E CO 
C. 1998701 IE 00 
C.20167008E 00 
C. 203870C6E 03 
C. 2C5870C3E 00 
0, 20787001E 00 
C. 20986998E 00 
C. 211869966 00 



0. 13757691E-C1 
0.14429174E-01 
0. 151 19277E-01 
0. 15 Q 27950E-01 
0. 1655513GE-01 
0. 1 7 3007325-01 
0. 18064659E-01 
0. 1 88468 10 E-0 1 
0. 19647073E-0 l 
0. 20465329E-01 
0.2130 1460E-0 1 
0.22 155326E-0 1 
0.23026783E-01 
0. 23915663E-01 
0 .249217926-01 
0.257449P2E-01 
0. 26665044E-01 
0.27641781 E-0 1 
0.286H990E-01 
0.29604476E-01 
0. 33610032E-01 
0. 31631447E-01 
C. 326685095-01 
0 . 33720996F-01 
0. 34798663F-01 
0. 35871338E-01 
0. 36966727E-01 
0.330806105-01 
0. 2 9 206754E-0 1 
0.403469175-01 
0.41500859E-01 
0.426633436-91 
0.43349129F-01 
0.45042980E-01 
0.46249662E-01 
0.47^63934E-01 
0. 48700567E-01 
0.49948323E-01 
0 .511999696-01 
0.52467272E-01 
0. 53745996E-01 
0. 550359126-01 
0.56336783E-01 
0. 57648379E-01 
0. 53970470E-01 
0. 60302628F-01 
0. 61645232E-01 
0. 62997460E-01 
0. 64359307E-01 
0.657305726-01 
0.6711 10 15E -01 
0. 63500519E-01 
0. 69898844E-01 
0 .713057526-01 
0.72721 1 24E-0 1 
0.74144721E-01 
0. 75576365E-0 1 
0. 770158176-01 
0. 78462899E-01 
0.7991 7431 E-0 1 
0. 81379235E-01 
0. 82948191E-01 
0.84324062E-01 
0.8 5 30672 7E-0 1 
0. 872960095-01 
0. 837917286-01 
0.902°3765E-01 
0.9180 2001E-01 
0.93316257E-01 
0. 94836414E-0 1 
0. 96362293E-01 
0. Q 795371 5E-01 
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C. 213 R 6993F 


00 


0.994 30561 E- 


-01 


0.21 5869 Q 0F 


00 


0. 13097271E 


00 


C.217P6998F 


00 


0. 10251999F 


00 


C. 21 986985? 


00 


G. 10407227F 


00 


C. 221 P6 5836 


00 


0. 10562944F 


00 


0. 22356980? 


00 


0. 10719138F 


00 


0. 22596977? 


00 


0. 10975791 E 


00 


0.22766975? 


00 


0. 11032397E 


00 


C. 22586972? 


00 


0. 11190450E 


CO 


0.23186970c 


00 


0. 11343432F 


00 


C. 233369675 


00 


G. 11506938E 


00 


0.23566965? 


00 


0. 11665654E 


00 


0.237869625 


00 


0. 1 1824870E 


00 


C.23986959E 


00 


0. 11984479E 


00 


0.241 869575 


00 


0. 12144470E 


00 


0.24336 554? 


00 


0. 12 30483 IE 


00 


0.2^596952? 


00 


0. 1 24 65 54 9 E 


00 


0.247863^9? 


00 


0 • 1262661 8F 


00 


C. 24986947? 


00 


0. 12799033E 


oc 


0.25186 c 445 


00 


0. 12949777E 


00 


0. 2530694 IE 


00 


0. 13111842? 


00 


0*25? q 6639E 


00 


0. 13274223E 


00 


C. 25 7 86 936? 


00 


0. 134369C8E 
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C. 71 245044E 00 
0.71419460E 00 
0.71 593875E 00 
0.7 1768290E 00 
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C. 93386066E 
C. 93=860636 
C. 93 786061 F 
C.93986058E 
C.941S6056E 
C. 94386053F 
C. 94RR6051E 
C. 94 7 860^RF 
0. 94 986 045 E 
C.95186043E 
0. Q53B6040E 
0.^5596038? 
0. 95786035E 
C.95986C32E 
C. 9618603 0E 
C. 96386027E 
r.965P6025c 
C.96786022E 
C. 9696602CF 
C.971P6C1 7E 
0. 97386014? 
C. -75860120 
0. 97766009 = 
C. 97986007 = 
0.931 96004F 
0.9938o0C2P 
C. 98595 9993 
C. 96785 9C6E 
0. 969 35994E 
0.991 65991 E 
C. 99335 989E 
C. 9<55P5986E 
C. 99735984E 
C.999B5981E 
C. 10013597E 
0. 1 003 85 95 E 
0. 10058594E 
0. 10078552E 
C. 10093591E 
0. 10113539E 
C. 101385 8 8E 
C. 10158587E 
C. 101785 85 E 
C. 10198584= 
C. 102185 82 E 
0. 10238581E 
C. 102585796 
C. 10278578E 
0. 10298576E 
0. 1031 8575E 
C. 10338573c 
0. 103585726 
0.103785716 
0. 103935896 
C. 10413566E 
0. 10438566 = 
C. 10458565E 
C. 1047 3563E 
C. 10498562E 
C. 105185o0E 
0. 1053P559E 
C. 10558558E 
0. 10578556E 
0. 10598555E 
0. 10618553c 
C. 10638552E 
C. 10658550E 
C. 106785^96 
0. 10658547E 
C. 10719546E 
0. 10738544F 
0. 10758543E 



00 


0. 719427C5E 


00 


00 


0.72117120= 


00 


00 


0. 72291535E 


OC 


00 


0.72465950= 


00 


00 


0. 7264C365E 


OC 


cc 


0. 72914780E 


00 


CO 


0.729P9196E 


00 


00 


0. 731t3611E 


00 


GO 


0. 73333026E 


00 


00 


0.7351 244 1 E 


OG 


00 


0. 7368685 6F 


00 


00 


0. 73861271F 


00 


00 


0. 74035686E 


00 


00 


0. 7421C101E 


00 


00 


0.7438^516= 


00 


00 


0.7^553932= 


00 


00 


0.74733347c 


00 


00 


0.7^907762= 


00 


00 


0.75082177E 


00 


00 


0.75256592E 


00 


oc 


0. 754310C7E 


00 


00 


0 . 75605422E 


oc 


00 


0.75779837E 


00 


00 


0.75954252E 


00 


00 


C. 78128o68E 


00 


no 


C.76303G83E 


00 


00 


0.76477498E 


00 


00 


0. 7665191 3E 


OC 


00 


0. 76326328E 


00 


00 


0. 7790C743F 


00 


CO 


0.77175156c 


00 


00 


0.77349573F 


00 


CO 


C. 77523988= 


00 


00 


0. 776QP404E 


oc 


01 


0.77872319E 


00 


01 


0. 73047234E 


00 


01 


0.73221 649E 


00 


Cl 


0. 73396064F 


00 


01 


0. 73570479E 


00 


01 


0. 73744R94E 


00 


01 


0.73919309E 


00 


01 


0.79093724F 


00 


01 


0.79268140E 


00 


01 


0. 79442555E 


00 


01 


0.79616970= 


00 


Cl 


0.79791 3 85F 


00 


Cl 


0. 79365800E 


00 


01 


0. 8Q140215E 


00 


01 


0. 80314630E 


00 


Cl 


0. 80489045E 


OC 


Cl 


0 . fi 066 3460E 


00 


Cl 


0.80R37876E 


00 


Cl 


0. 81012291E 


00 


01 


0. 811R6706E 


00 


01 


0.81361121E 


00 


01 


0. 81535536F 


00 


01 


0. P1709951E 


00 


01 


0.91384366E 


00 


01 


0.820587816 


00 


01 


0. 822331 C 6E 


00 


01 


0.82407612E 


00 


01 


0. 8258202 7E 


00 


01 


0 . 82 756442E 


00 


01 


0. 82930857F 


00 


01 


0. 83105272E 


00 


Cl 


0. 8327 Q 687E 


00 


01 


0.63454102= 


00 


01 


0.83628517= 


00 


01 


0. 83802933E 


00 


01 


0.83977348E 


00 


01 


0.84151763E 


00 


01 


0. 84326178E 


00 
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C.107785426 Cl 
0.10798540? 01 
0.108185396 01 
C, 10 839537E 01 
C. 108565366 01 
0. 10878534E 01 
C. 10698533E 01 
0.109185316 01 
0.10938530c 01 
0.109585296 01 
0.109785276 01 
C. 109985266 01 
0. 110185 24 E 01 
0. 11038523E 01 
0.110585216 01 
C. 110785206 01 
0.110985186 01 
C. 111185175 01 
0 . 11 1385156 01 
0.111585146 01 
0.111785136 01 
0.111985115 01 
0.112185106 01 
0.112385086 01 
0.112585076 01 
0.112785056 Cl 
0.112965046 01 
0.113185026 01 
0.113385016 01 
0.11 358500E 01 
0.113784986 01 
C. 11398497E 01 
0.114184956 01 
0.11438494? 01 
0.114584926 01 
0.114784916 01 
C. 114984896 01 
0.115184886 01 
0.11 538486E 01 
0.1155 3485E Cl 
C. 115784846 01 
0.115984826 01 
C. 116184816 01 
C. 116384796 01 
0.116584786 01 
0. 1167 3476E Cl 
C.11698475E 01 
0.1171 84 73E 01 
0.11 738472E 01 
C.11758471E 01 
0.11778469E 01 
C.11798468E 01 
0. 11818466E 01 
0. 1 1 838465F 01 
0. 1 1 85S463F 01 
0.118784625 01 
0. 11896460E 01 
0 • 11916459E 01 
C.11938457E 01 
0. 11958456E 01 
0.11978455E 01 
0.119984536 01 
0.120184525 01 
0. 120384506 01 
0. 12C584496 01 
0. 12078447E 01 
0. 12098 44 6E 01 

0.121184445 01 
0.121384436 01 
0.12 158442E 01 
0.121 78440E 01 
C.12198439E 01 



0.845005876 00 
0.846 74984E 00 
0.848493586 00 
0.850236956 00 
0.851979856 00 
0.85372216E 00 
0.8 5 546 36 8E 00 
0.857 2042 6E 0C 
0.858943766 00 
0. 86C68207E 00 
C.86241895E 00 
0.864154226 00 
0. 8653 8 77C 1 - 00 
0.867619226 00 
0.869348476 00 
0.87107527E 00 
0.872799286 00 
0.874520246 00 
0.876237876 00 
0.377951866 00 
0.879661866 00 
0.88136750E 00 
0.833068446 00 
0.884764315 OC 
0.93645470E 00 
0.838139196 00 
0.8398 1 742E 00 
0.89143891E 00 
0.993153256 OC 
0. 394810026 00 
0.896453756 00 
0.89309901E 00 
C.89973032E 00 
0.90 135223E 00 
0.90296429E 00 
0. 99456605E 00 
0 • 906 1 5702E 00 
0.90773672E 00 
0. 90930480E 00 
0.91086072E 00 
0.9 1240412E 00 
0.91393453E 00 
0.91545159E 00 
0.9 169548 2E 00 
0.91844386E OC 
0.91O91836E 00 
0.921377906 00 
0.922822185 00 
0 • 92 425084E 00 
0 . 92566353F OC 
0.92 706001E 00 
0. 92343997F 00 
0. 9298031 3E 00 
0.931 149 1 8E 00 
0 • 932477°5E 00 
0.933 7 892 5E 00 
0. 93508285E 00 
0.936353636 00 
0.9 376 1635E OC 
0.93385589E 00 
0. 94007713E 00 
0.94127995E 00 
0 • 94246429E 00 
0. 94363010E 00 
0. 94477731E 00 
0 • 94590592F 00 
0. 94701588E 00 
0.34810718E 00 
0.94917989E 00 
0.950233946 00 
0.95126945E 00 
0.95223642E 00 
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C.12218437E Cl 
0. 1223S436E Cl 
0. 1225 e 434 E 01 
0.122784335 01 
0.122984315 01 
0. 123184305 01 
0.123384295 Cl 
C. 12 2 5 342 7 6 Cl 
0.123784265 01 
0. 123984245 01 
C.12M8423E 01 
C.12438421F 01 
C.1245842CF 01 
0. 1247941 9E 01 
0.12^984175 01 
C. 125184155 01 
C. 125394145 01 
0.12558413E Cl 
0.12579411E 01 
C.12599410E Cl 
0 . 12613409E Cl 
C.12639407F Cl 
C. 126584055 01 
0. 1267940^5 01 
C. 126334025 01 
0 . 127184C1E Cl 
C. 12 7384C0E Cl 
0. 127583965 01 
C. 127793975 01 
C. 12748395c 01 
0. 128183 94 f 01 
0 . 12638392E Cl 
0 • 12853391E Cl 
C.12878389E 01 
C.1263«383F Cl 
C.12918386E Cl 
C. 1293S335E Cl 
0. 1235b334 c 01 
C • 12 C 76382E 01 
C. 12998381E Cl 
C.13C18379E 01 
C . 1 3038378E 01 
C. 13C59376F 01 
0. 1 3078375F 01 
C. 1 30993 73 E 01 
0.131193 72 E 01 
C.13138371E 01 
0. 13153369F 01 
C. 131783635 01 
0 • 131933665 01 
0. 13218365E 01 
C. 13238^635 01 
0. 13259362E Cl 
C. 132783605 01 
f. 132983 59F 01 
C.1331»357E Cl 
C. 13338356E Cl 
i" 1 . 1 33283555 01 
0. 1337 9353E Cl 
C.13388352E 01 
C.13418350E 01 
0.13^383495 01 
0. 134 5834 7E 01 
0. 134-73346E 01 
0.134983445 01 
C.13518343E 01 
0. 1 3530342E 01 
0. 1 3558340E 01 
0. 13578339E 01 
C. 1 3598337E 01 
0.13618336F Cl 
C. 1 3638334E 01 



0.95328496F 00 
0.954265245 00 
0.95522726E 90 
0.95617U5E 00 
0. 95709705E 00 
0. 95800507E 00 
0.95.9895395 00 
0.959768125 00 
0. 96062338E 00 
0.96146142E 00 
0.962282425 00 
0. 96308655E 00 
0.96387398E 00 
0.96+644 c 15 00 
0.96539956F 00 
0.966138185 00 
0. 96686095E 00 
0. 967568045 00 
0.968259695 00 
0.969936155 00 
0.9695976^5 00 
0.9702^415 00 
0.97097669E 00 
0.971494735 00 
0.972098775 00 
0.97263897E 00 
0.9?32 A 565E 00 
0.973829035 00 
0. 97437936E 00 
0. 97491698E 0C 
0. 9754^1 75E 00 
C.97595429E 00 
0.97 6^54685 00 
0.9769431 9E 00 
C.97742003E 00 
0.977895485 00 
0.97833973F 00 
0.9797830 IE 00 
0.97921556E 00 
0.97963762F 00 
0 . 9 9004 9435 00 
C. 990451175 00 
0. 990943075 00 
0.93122537F 00 
0. Q 8 1598265 00 
0.991961915 00 
C.98231655E OC 
0.93266238E 00 
0.93299 9 56E 00 
0.°3332834E 00 
0.93 3648°0E 00 
0.933961^0E 00 
0. 9342661 C5 00 
0.934563175 00 
0.994952735 00 
0.935 13502E 00 
0.935410215 00 
0.985678495 00 
0 . 93593998E 00 
0.936194855 00 
0.936443225 00 
0.936685285 00 
0.93692119E 00 
0.937151 09E 00 
0.987375085 00 
0.98759335E 00 
0.937306025 00 
0.983013275 00 
0.93821527E 00 
0. 93841208 E 00 
0.939603955 00 
0. 93379093E 00 
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C. 1365 8 3336 Cl 
0.13673331E Cl 
0. 1 35 9333CE Cl 
C. 137193286 01 
C. 137383276 Cl 
C. 1 3758326E 01 
C.13778324F 01 
C.137CR323E 01 
C. 138183216 01 
0. 1 383832CE 01 
C.13858318E 01 
0 . 136783176 01 
0.138983156 01 
C.13918314F 01 
0.139383136 01 
C.13958311C 01 
0. 1 3978 3 1 C c 01 
0. 13996308F 01 
C. 140133 07E 01 
0*140383056 01 
C.140583C4F 01 
0. 140733026 01 
C • 1409 8 3016 01 
C. 141182996 01 
0. 141382985 01 
0.141532976 01 
C. 141782955 01 
C . 1 4 1 <5 92 946 01 
0.142182926 01 
0.142382916 01 
0. 14253289E 01 
0 • 14273288E Cl 
0.142982866 Cl 
C • 1 43 1 82 85 F Cl 
0. 1 4333284F Cl 
C. 143502825 Cl 
C. 14378281E Cl 
C.14398279E Cl 
C. 144182786 01 
0.144382765 01 
C.14453275F 01 
C. 14478273E 01 
C. 144932726 01 
0, 145 1 8270E Cl 
C. 145382696 01 
C. 14 5582686 01 
0. 1457826 6F 01 
0. 145982656 01 
0.146182636 01 
0.14638262F 01 
0.146582606 01 
0.146782596 01 
0.146982575 01 
0.147182566 01 
0.147382556 01 
C. 147582536 Cl 
0. 147732526 01 
C. 147982506 Cl 
C. 148182496 01 
C. 148382476 01 
0. 148582466 01 
0.148782446 01 
0.148982436 01 
0.14^132426 01 
0.149382405 01 
0.149582396 01 
C. 149782376 Cl 
C. 149982366 01 
C. 150132346 01 
0. 15C38233E 01 
G.15C58231E 01 
0.15073230F 01 



0.93897314E 00 
0.93915070E 00 
0.93932374E 00 
0.939492366 00 
0.939656696 00 
0.9898 1^786 00 
0.939972776 00 
0.990124766 00 
0.990272826 00 
0.990417066 00 
0.990557616 00 
0.990594526 00 
0.99082792F 00 
0.990957926 00 
0.9910 q 464E 00 
0.991203146 OC 
0.991323466 00 
0. 99144578F 00 
0. 991560 16E 00 
0.95167162F 00 
G.99 173028F 00 
0.991866146 00 
0.99196 C 31P 00 
0.992089876 00 
0.9921 9786E 00 
0.99228334E 00 
0.992376396 00 
0.992467056 00 
0.992555386 00 
0.992641456 OC 
0.992725316 00 
0.992807036 00 
0.992886666 00 
0.992964276 00 
0.993039976 00 
0.99311376F 00 
0.993135706 00 
0.993255856 00 
0.993324286 00 
0.993390966 00 
0. 99345601 E 00 
0.99351 94 36 00 
0.993581246 00 
0.9936415G6 00 
0.993700276 00 
0.993757556 OC 
0.99381340E OC 
0 . 9 93 867 8 1 E 00 
0.99392086E 00 
0.993972546 00 
0.994022916 00 
0.994072026 00 
0.994119886 00 
0.994166556 OC 
0.99421203E 00 
0.994256386 00 
0.99429965E OC 
0.994341856 00 
0.994382986 00 
0.99442309E 00 
0 . 9944622 5F 00 
0.994500466 00 
0. 99453777F 00 
0. 99457413E 00 
0.994609596 OC 
0.99464422E 00 
0.994677966 00 
0. 9947108 6F 00 
0.994742996 00 
0.99477^285 00 
0.994804306 00 
0.994834546 00 
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C. 1 5C98228E 
0. 15118227E 
C. 15138226E 
0.151*82246 
0.151732236 
C. 151 C 8221F 
r. 15218220F 
C. 1523821PF 
C. 1525S217E 
0. 1*2782155 
0.152982145 
0. 153182136 
C. 153382115 
0. 15359210E 
0. 1 5378208E 
C. 153902076 
0. 154182C5E 
0. 1?^382C4E 
0. 154532026 
C. 15478201E 
0. 1549°199F 
C. 155181986 
C. 1 C 538197E 
C. 15558195 E 
C. 155781945 
0. 15 59 8192 E 
0. 18618191E 
0.15638189E 
C. 156581886 
C. 156791865 
0. 15698185E 
0. 1571S134E 
0. 15738182E 
0. 157581816 
O. 15776179E 
0. 15798178E 
0.156181755 
0.158381755 
0. 15856173E 
0.158781725 
C. 15898170E 
C. 15918169F 
C. 15933168c 
0. 15 C 58166E 
0. 15978165F 
0. 15998163E 
C. 16013162E 
G. 16C38160E 
C. 15C58159F 
0. 1607615 7 E 
0. 16098156E 
C . 1 6 1 1 8 1 5 5 E 
C. 1612ei?3£ 
C. 16163152E 
0. 16173150 c 
0. 16138149E 
0.162181475 
C. 16233146F 
C. 1625814AF 
0. 162781435 
0. 16298141E 
0. 163181405 
r. 163381395 
C. 163531375 
0. 16378136E 
0.163991345 
C. 164 1 31 3 3E 
0. 164381315 
0. 1645813C5 
0. 164791285 
C. 164981275 
0. 165131265 



0.99486357E 00 
0.99489188E 00 
0.99491Q46F 00 
0.994946366 00 
0.994972595 00 
G. 994998165 00 
0.995023075 00 
0.995047395 00 
C.99507105E 00 
0.9 95094 1 8E 00 
0.99511671F 00 
0.995138705 00 
0.9951 60 16E 00 
0.99519108E OC 
0.99520153E 00 
0.99522144F 00 
0.995240875 00 
0.995259826 00 
0.99527836F DO 
0.995 2 C 642E 00 
0.995 314065 00 
0.99533129E 00 
0.995 39 810E 00 
0.99536449E 00 
0*99 5 3809 6E 00 
0.995396C2E 00 
0.995911165 00 
0.995925946 00 
0.995440365 00 
0.93545943F 00 
0.99546*146 00 

C.9954814CP 00 

0 . 9 954 9448 E 00 
0.99 55071 2 E 00 
0.995519465 00 
0.995531505 00 
0.99554318E 00 
0.99555457E OC 
0.99556565E OC 
C.9955765CE 00 
0.995597C5E 00 
0.99559730E 00 
0.9956C732E 00 
0.99561709E 00 
0.995626636 00 
0.99563593E 00 
0.935 64499E 00 
0.99565381E 00 
0 • 9 °5 662 3 9E 00 
0.9 3567080E 00 
0.935678°6E 00 
0.99568695E 00 
0.99569470E 00 
0. 9957022 7F 00 
0.995709606 00 
0.99571675E 00 
0.99572372E 00 
0.99573052E 00 
0. C 9573714E 00 
0.99574357E 00 
0.995743835 00 
0.99575 59 IE 00 
0.99576181F 00 
0.9 9576753F 00 
0.995773086 00 
0.33577844E 00 
0.99579369* 00 
0.99578*755 OC 
0.99579370E 00 
0.99579847E 00 
0.99580312 c 00 
0.99580759F 00 



01 

01 

01 

01 

Cl 

01 

01 

01 

01 

Cl 

01 

01 

Cl 

Cl 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

Cl 

Cl 

Cl 

01 

01 

01 

Cl 

Cl 

Cl 

Cl 

01 

01 

01 

01 

Cl 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

Cl 

01 

01 

01 

01 

01 

01 

Cl 

01 

01 

Cl 

01 

Cl 

Cl 

01 

01 

01 

01 

01 

01 
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C . 165331246 
0. 164=18123?. 
C. 165781215 
0. 16598120E 
C. 16613118? 
0.16633117? 
0.166581158 
0. 166781 14E 
0.16698112E 
n. 16718111E 
0. 1 6 7 3 3 1 1 Of 
0. 167581C8E 
0. 167781C75 
0. 167661C5E 
0. 168131 0=6 
C. 16P331C26 
C. 16853101? 
C. 168780996 
C . 16893056E 
n. 16518097F 
C. 169380950 
C. 165590946 
C. 169730026 
0. 169930916 
C. 1 70100896 
C. 17C30O83F 
0. 17C56086E 
C. 17C780856 
C. 170980836 
C. 17118C82F 
C. 171380616 
0.171530796 
f .171780786 
C. 171980766 
0. 1 72 1 8075 6 
C. 172380736 
0.17258072E 
0. 172730706 
0. 17298069E 
C. 173130686 
0.173380666 
0.173530656 
C. 17373063E 
C. 173980626 
0. 1 74180606 
0. 1743 3 0596 
0.17453057 F 
0. 174780566 
0.174=30556 
0.175180536 
0. 17538052E 
0. 175580506 
0.175780496 
C. 175930^76 
0.176180466 
0. 176380446 
0.176580436 
0.176730416 
C. 176980406 
0. 17718039E 
0. 1773B037E 
0. 177580366 
C. 177730346 
0. 177QS0336 
0. 178130316 



01 

01 

01 

01 

Cl 

Cl 

Cl 

Cl 

01 

01 



01 

01 

01 

Cl 

Cl 

01 

Cl 

01 

01 

Cl 

Cl 

Cl 

01 

Cl 

01 

01 

Cl 

Cl 

Cl 

01 

01 

01 

Cl 

01 

Cl 

Cl 

01 

Cl 

01 

01 

01 

Cl 

01 

01 

Cl 

01 

01 

Cl 

01 

Cl 

01 

01 

01 

Cl 

01 

01 

01 

01 

01 

01 

01 

01 

Cl 



0.995311946 00 
0. 9958161 7E 00 
0.995820286 00 
0.995824286 00 
0.995828096 00 
0.995331796 00 
0.99583536E 00 
0.995839886 00 
0.9958422*6 00 
O.99504556E 00 
0.99584 372E 00 
0.995851766 00 
0.99535^746 00 
0.995857606 00 
0.995360406 00 
0.99536308 F 00 
0.995465646 00 
0.995885156 00 
0.99587053F 00 
0.995072866 00 
0.995 37512 p 00 
0.995877276 00 
0.995379356 00 
0.995331326 00 
0.995383236 00 
0.995085076 00 
0.995986806 00 
0.995898476 00 
0.995890086 00 
0.995391576 00 
0.995893006 00 
0.995894376 00 
n . no 



0.995866946 OC 
0.995893136 00 
0.95569=206 OC 
0.995900216 OC 
0.995901176 00 
0.99590206E 00 
0.9959C290E 00 
0.99590367E 00 
0.995904396 00 
0.995905046 00 
0.995905646 00 
0.995906176 00 
0.9=5906656 OC 
0 • 9 95 90707E 00 
0.995907496 00 
0.595907646 00 
0.99590814= 00 
0.99590838E 00 
0.995904566 00 
0.995908746 00 
0.995906866 00 
0.995908926 00 
0. 995900=86 00 
0 . 99590898E OC 
0. 99590898E 00 
0. 995903986 00 
0.99590898E OC 
0.995908986 00 
0.995908986 00 
0.995908986 00 
0. 995908986 00 
0.995908986 00 
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APPENDIX II 



The tabulated sample data for the non-dete rmmistic case, Table II, 
was obtained from a modified form of equation (19). Using the 
modified equation. 



d) = 0. 175 + i 

Y Y(y) X(x) 

where X(x) is N(130;10) 

Y(y) is N(65;12) 

values of 0 and sensitivity were calculated for chosen values of x and 
y. Assuming these to be measured values, the data was then used to 
obtain the cumulative distribution of 0. 

Since the data indicates 0 is more sensitive to changes in y than 
to changes in x, the problem is worked with x as a constant and y as 
a single variable as discussed in section III. Aitkin-Hermite inter- 
polation [Hildebrand 1956] is used to estimate the values of y between 
data points. This interpolation scheme utilizes the sensitivity data 
and obtains a more accurate estimate than would linear interpolation. 
The calculation of the cumulative distribution of 0 is then made 
using equation (35) in discrete form. As in the deterministic case, 
a Ax of 0. 002 was used for rectangular integration. 
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1 PP END I X IT 



PROGRAM ONE 



C CALCULATION OF THE CUMULATIVE DISTRIBUTION OF AN 

OUTPUT, PHI , FROM MEASURED VALUES OF THE PARAMETERS 
X, Y, AND THE SENSITIVITY OF PHI WITH RtSPECT TO X 
AL m Y. THE PROBABILITY DISTRIBUTIONS OF X AND Y ARE 
KNOWN, N(lBOSlO) AND N(65;12> RESPECTIVELY, AND THE 
DATA POINTS HAVE BEEN SELECTED IN DX INTERVALS TO 
FACIIITATF THE CALCULATIONS. 

D I Mf NS I ON ARGI ( 5,5) , VAL( 10 > , ARO( 5 ) ,XC (5 ) ,F ( 5, 2 ) , Z ( 5) , 
FYT( 25 ) , PXC( 5) , XI (200) ,PX (200 ) , YT( 25) 

C READ THE MEASURED VALUES OF X AMD THEIR CORRESPONDING 

C PFOPAPT LITIFS, DENOTED BY XC AND PXC. 

READ (5, 102) ( XT ( I ) , P XC ( I ) ,1 = 1,5) 

WRITE (6, 1C?) 

10? FCRMAT(//,T21, 'XC' ,T47» 'PXC' ,/) 

WF I Tfc ( 6 , 1 Cfe ) ( XC ( I) , PXC ( I ) , 1 = 1 , 5 ) 

106 FORMAT ( 2F26. 5 ) 

101 FCRMAT(F10.5) 

REAP THE MEASURED VALUES OF PHI , DENOTED BY ARGI. 
THERE IS A 15 0 1 N T PHI CORRESPONDING TO EACH SET OF 
X - Y POINTS. 

»t ACM 5 , 101 ) ( ( ARGI ( I , J) , J = 1 ,5 ) , 1 = 1,5) 

C READ THE MEASURED VALUES OF Y AND THE CORRESPONDING 

SENSITIVITY OR PHI WITH RESPECT TO Y FOR EACH OF 
THE S c VALUES, DENOTED BY F(I,J). 

PC AIM f , 102) ( ( P( I , J) ,J = 1 ,2) , 1=1 ,5) 

102 FCPMAT(2F10.5) 

F E AC THE KNOWN CUMULATIVE DISTRIBUTION OF Y, DEMOTED 
BY YT AND PYT. 

READ ( 5 , 1C2 ) ( YT( I ) ,PYT( I ) , 1 = 1 ,25 ) 

XI ( 1 ) =0. 5708710 
DC 20* K=1 ,168 
P T=C. C 

C X DENOTES THE VALUE OF PHI FOR WHICH THE CUMULATIVE 

C DISTRIBUTION IS BEING CALCULATED. 

X=X1 ( K) 

DC 2C2 1=1,5 
DC 20 1 J = 1 » 5 

Z(J) DENOTES THE ARRAY CF ARGI VALUES CORRESPONDING 
TC A PARTICULAR VALUE OF XC. 

201 7 ( J)=APG1 ( I , J) 

DX=1 C. 0 
F R S= 0 . CC 1 
Nr IM=5 

C SUBROUTINE ATSM ORDERS THF ARRAY Z(J) SUCH THAT 

ABS(APG( I )-Xl.GE.ABS(ARG(J)-X) IF I.GF:J. THIS IS 
REQUIRED FOR SUBROUTINE ONE. 

CALL ATSM(X,Z,F,5,2,ARG,VAL,NCIM) 

C SLBRGUTINF ONE I NTERPCLATES c UNCTION VALUE Y FOP A 

GIVEN ARGUMENT VALUE X USING THE TABLt ( A R G , V A L ) OF 
ARGUMENT, FUNCTION, AND DERIVATIVE VALUES FRCM THE 
C SUBROUTINE ATSM. THIS SUBROUTINE USES THE AITKEN- 

HtFMITE I R TERPGL AT I ON. 
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CA LL QMRI X,ARG f VAL,Y,NDIM,EPS , I £R ) 

C SUBROUTINE TWO USES LINEAR I NTE P POL AT I CN TC OBTAIN 

P ( Y ) FOR The VALUE OF Y OBTAINED IN SUBROUTINE ONE . 

CALL TW0(Y,PY,YT,PYT) 

PT=PT+< I-FY) *PXC ( I ) 

P X ( K )=pt*qx 
332 FCR M AT(2E3C.8) 

2 0? Cf NT INI 'F 

X1(K+1) =X1 ( K ) + 0. 0020000 
204 CCNTINUF 

WF T TF < i, 107) 

1C 7 FCPNAK 1HI ,T20 U( I )• ,T42 , 'CUMULATIVE DISTRIBUTION' , / / 
WFITF(t,3C?) ( X l ( K ) * PX ( K ) , K= 1 , 1 6 9 ) 

STOP 

ENT 



SI P P OUT INF ATS*MX,Z,F, I ROW, IF CL , ARC , v AL , NC I M ) 

D IMF NS I CL Z ( 1 ) , F ( 1 ) ,ARG( 1 ) , VAL ( 1 ) 

C CASE I c PW = 1 IS CHECKED CUT 

IF< IFOV -1 ) 2 3 . 21 ,1 

1 f i=nt i m 

C IF N IS C c E A Tc p ThAN IRDW, N IS SET EQUAL TC IRCW. 

I F< N-IkPW ) 3 t 3, 2 

2 N= I P CV 

C CASE IRCW .GE. 2 

C SFAPCHING FOR SUBSCRIPT J SUCH THAT Z(J) IS NEXT TO X. 

3 I F ( Z ( I p PW ) - Z ( I )) 5,4,4 

4 J = I R f W 
1 = 1 

GP Tr 6 

5 I=IPPW 
J =1 

6 K — ( J + I ) /? 

IF (X-Z ( K ) ) 7 , 7, 8 

7 J = K 

GC T c 
a I =K 

7 IF ( I A p S ( J-I ) -1 ) 10, 10,5 
10 IF(APS( 7 <J)-X)-ABS(Z(n-X))l?, 12 ,Ll 
1 1 J = I 

C TABLE 5 El EC T I ON 

12 K = J 

J L =C 

jr<=o 

DC 20 I = 1 , E 
AFG(I) =Z ( K ) 

I F < ICCL-l) 14 , 14,13 

1 3 VAL( 2 *I- 1 )=F<K) 

kk=k+irpw 

VA L ( 2* 1 )=F(KK) 

GC TP 1 ? 

14 VALI I) =F ( K ) 

15 JJP-J+JR 

IF( JJR-IPPW) 16 , 18, 18 

16 J JL= J- JL 

I F ( JJL -1 lie, 19, 17 

17 IF ( ARS( Z ( JJP+1 ) -X)-ARS< Z( JJL-1 )-X ) >19,19,18 

18 JL =J L 4 1 
K= J- JL 
GC TO 2 C 

19 JF =JP + 1 
K = J + JR 

20 CONTINUE 
RF TURN 
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C CASE IFCW=1 

21 AFG( 1 )=7 ( 1 ) 

VAL( 1 ) = F ( 1 ) 

I F ( 1 C OL- 2 ) 2 3 » 22 » 2 3 

22 VAL< 2 )=F (2) 

23 P F T URM 
END 



SURF OUT I NF ONE ( X f ARG , VAL , Y , NO I M , EPS , I ER ) 

D I M FNSI ON A RG ( 1 ) ,VAL( 1 ) 

IFP = 2 

H2=X-ARG( 1 ) 

I F ( ND I M- 1 ) 2 » 1 1 3 

1 Y=VAL( 1)+VAL(2)*H2 

2 RETURN 
C 

C VECTCK APG HAS MORE THAN 1 ELEMENT • 

C THE FIRST STEP PREPARES VECTOR VAL SUCH THAT A I TK EN 

SCHEME CAN RE USED. 

3 1 = 1 

PC 5 J = 2 »NDI M 
H 1 =H 2 

H2=X-ARG( J ) 

Y= VAL ( I ) 

VAL ( I )=Y+VAL ( 1+1 ) *H1 
H=H 1-H2 

I F < H ) ^ , 1 3 , A 

4 VAL( 1+1 )=Y+< VAL( I+2)-Y)*Hl/H 

5 1 = 1+2 

VAl ( I ) = VAL ( I ) + VAL ( I + i )*H2 
C END OF FIRST STEP 

C PREPARE AITKEN SCHEMP 

PE LT2=0. 

I I ND = 1 - 1 

C START A I TKEN-LCCP 

DC c I = 1 , 1 END 
PELT1=PELT2 
Y=VAL ( 1 ) 

M=< I+3J/2 
H 1=AP G ( M ) 

DC 6 J = 1,I 
K= I + 1- J 
L=(K+l)/2 
H=ARG ( L ) -H 1 
IF(H)fc,l* t 6 

6 VAL (K ) = < VAL ( K ) *< X-Hl )-VAL (K + l )* < X-ARG(L )) > /H 
DELT2=ARS(Y-VAL( 1) ) 

I F ( DE L T2-E FS ) 11 ,11 ,7 

7 ir(I-5)9,8,3 

3 IF(rtLT2-Df LT1)R,12,12 
9 CONTINUE 

C END OF AITKEN-LCOP 

10 Y=VAL( 1 ) 

R F T UF N 

THERE IS SUFFICIENT ACCURACY WITHIN 2*NCIM-2 ITERATION 
S T E F $ 

11 I FP=C 
GOTO 10 

C TEST VALUE DELT2 STARTS OSCILLATING 

12 I CR=1 
RETURN 

C THERE APF TWO IDENTICAL ARGUMENT VALUES IN VECTOR A°G 

13 Y = VAL( 1 ) 

14 I ER=3 
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KE TUP N 
FND 



S( QF-PUT I NF TWP( Y ,PY,YT,PYT) 

DIMENSION V T ( 2 5 » ,PYT(25) 

IF ( Y.LT.YTE 1 ) ) GO TO 301 
J = 2 

110 IFIY-YTIJI) 200,201,202 
202 J=J+1 

IFIJ.LF.25) GO TO 110 
P Y= 1 « 0 
GC t c 300 
201 PY=PYT(J) 

GC TC 300 

23C PY = FYT(J-l) + (PYT(J>-PYT(J-in/(YT(J)-YT(J-in*(Y-YT(J- 

GC TC 3 CO 
30 l P Y =C . 0 
TOC F f TURN 
END 
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PHI 

C.57C871CCE 
0. 572S7097F 
0. 57497094E 
C.57687092E 

O. 57887C89E 

P. 53037087E 

C. 582 87C8A? 

0. 58487082E 
0. 58687C79E 
0. 58387C76E 
0. 5908707^E 
0.59287071E 
0. 594P7C69E 
0. 59667066E 
C. 598870648 
0* 60087061E 
C.60287053E 
C • 6 04 8 7 0 5 6 E 
0. 60687053E 
C. 6C387051E 
C.61C87048E 
0. 61287C45E 
C. &1487043E 
0. 61687040E 
C. 61 98703SE 
C.62C87035E 
0. 62287033E 
0. 62487030E 
C.62687027F 
0. 62887025E 
0.63C87022E 
0. 63297 02CE 
0. 6348701 7E 
0. 6368701 5E 
0. 63887012E 
0. 640870C9E 
C.64287007E 
0. 64467004E 
0. 64687002E 
C. 64886959E 
C. 65C86997E 
0. 65286994E 
0. 65486 9 91E 
C. 65686989E 
C. 65 886986E 
C. 66036984E 
0. 66286981E 
C.66486979E 
0. 66686976E 
0. 66986973E 
0. 67086971E 
0. 672 86968E 
C.67486966E 
0. 67686 96 3E 
0. 67 896961 E 
0.68C86958E 
0. 68286955E 
0. 68486953E 
0. 6868695CE 
C. 68886 948E 
0.69C86945E 
0. 69286942E 
C. 6 9486940E 
C.69686937E 
0. 69886935E 
0. 70086 932E 
0.70286930F 
0. 70486927E 
0. 70686924E 



CUMULATIVE 





C I STP I BUT ION 


00 


0 • 62946975E-01 


00 


0.681461 1 OF-O 1 


00 


0.74001014E-01 


00 


0.797 9 868 5E-0 1 


00 


0. 85879432E-01 


00 


0.92248380F-01 


00 


0.98956390E-01 


00 


0.10555303E 00 


00 


0.11299674F 00 


00 


0. 12032640E 00 


00 


0.1278770CE 00 


00 


0. 13566995E OC 


00 


0.14363253E OC 


00 


0.15177727? 00 


00 


0.16019422'= 00 


00 


0. 16903991E 00 


00 


0. 17780799E 00 


00 


0. 19684916E OC 


00 


0.1 9600827E 00 


00 


0.20516866E 00 


00 


0.21465921E 00 


oc 


0 • 224080 1 5E OC 


00 


0.23413759E 00 


00 


0 • 2440890 IE 00 


00 


0.25411922E OC 


oc 


0 • 264 3328 9F 00 


CO 


0.27449834E OC 


00 


0. 2346822 7E 00 


00 


0.29506922E OC 


00 


0.30539531E 00 


00 


0.3161C376E 00 


00 


0.326719466 00 


00 


0. 33738047E 00 


oc 


0.349 10 561E 00 


oc 


0. 35374993E 00 


00 


0 . 369 3 3064E 00 


00 


0.37996292 c 00 


CO 


0 . 39054 102E 00 


00 


0.40114039E 00 


00 


0.411 77332F 00 


00 


0.42 2 294 1 4E 00 


00 


0 • 43 2 8643 IE 00 


oc 


0.443 32552E 00 


00 


0.45 37053 7E 00 


00 


0 . 46397537E 00 


00 


0.474 1 5489E 00 


00 


0 .484 2432 7E OC 


00 


0.49421632E 00 


00 


0.50409496E 00 


00 


0.51387882E 00 


00 


0 . 52 3 56750E 00 


00 


0.53316504? 00 


00 


0.54267335E OC 


oc 


0 • 55 209064E OC 


00 


0.56139952E 00 


00 


0. 57061529E 00 


00 


0. 57962447E 00 


00 


0 • 58854502E 00 


oc 


0.59 722644E 00 


oc 


0. 6057699 9E 00 


oc 


0.61422157E 00 


00 


0.6224 82 12E 00 


00 


0.63067019F 00 


00 


0.63879596E 00 


00 


0. 646801 59E 00 


00 


0.65472543E OC 


00 


0 • 66255569E 00 


00 


0. 6702089 3F 00 


00 


0.6776 9 872E 00 
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c. 


, 70B86922F 


00 


c. 


, 71086515 = 


00 


c. 


, 71 28691 7E 


00 


c. 


, 7143691*8 


00 


c. 


,71686912= 


00 


c. 


, 71 8869C9E 


00 


c. 


, 72C869C6E 


CO 


0. 


, 722869042 


00 


0. 


, 72486901 = 


OC 


0. 


726 6 685 9 E 


00 


0. 


, 728868 96E 


00 


c. 


,7303689*E 


00 


c. 


73286=91= 


00 


0. 


, 73*96888E 


00 


G. 


, 7 3686 8 86E 


CO 


c . 


73885393= 


00 


0. 


, 7A036 361E 


00 


c. 


, 742 8 53 78 E 


00 


6 1 


,744=6=76= 


00 


c. 


, 746=6873= 


00 


c. 


, 7488637CE 


00 


c. 


, 7503686RE 


OC 


0. 


, 75286 R555 


00 


c. 


, 75*863638 


00 


c. 


. 7568686CE 


CO 


c. 


, 7588685RE 


OC 


0. 


, 76086855E 


00 


0. 


, 76286 952E 


CO 


c. 


,76466350= 


00 


c. 


, 766S6B47E 


00 


c. 


, 76886845E 


OC 


c, 


.77086342= 


00 


c. 


, 772R6339E 


00 


c . 


,77436837= 


00 


c. 


,77686=34= 


00 


0, 


, 77886332E 


00 


c. 


, 7 8086 32 9E 


00 


e, 


,79286 3 27E 


00 


0. 


,7 8486 324F 


00 


c. 


, 73636821E 


00 


c . 


, 78886315E 


CO 


c. 


, 79 086 3 1 6E 


00 


c, 


, 792R6314F 


00 


0, 


, 79486311E 


00 


c. 


, 79686309F 


00 


0. 


, 793 = 6 306E 


00 


c. 


, R0C863C3E 


00 


0, 


, R0286B01F 


00 


c. 


, 304 86 7 C 8E 


00 


C. 


, 80685756= 


00 


0. 


, 30886793F 


00 


c. 


, 81066791= 


CO 


0. 


, 512667 88 E 


00 


c. 


, 81 486785E 


00 


0. 


,816867838 


00 


C i 


, 31 88678CE 


C 0 


0. 


, 82 0 86 778 E 


00 


c, 


,3226 5775E 


CO 


G, 


, B2A06773E 


OC 


G. 


, 82686770= 


00 


c, 


, 82=66767= 


00 


c. 


, 83C86765E 


00 


0. 


, 83286762E 


00 


c, 


, 83486760E 


00 


c. 


, 83686757E 


00 


c. 


, 33886755E 


00 


c, 


, 84C86752E 


OC 


c. 


, 84266749E 


00 


0. 


, 844S6747E 


OC 


c. 


, 846 86 744 E 


OC 


0. 


, 84886742E 


CO 


G« 


, R5086739E 


00 



0.634999C4E 00 
0.69199145= 00 
0. 69892645E 0C 
0.7056S224F 00 
0.71225584E 00 
0.71380460E 00 
0.725 27 1705 00 
0.73166728= 00 
0 . 7 3 30 1 1 60F 00 
0 . 7AA25936E 00 
0.75045109E 00 
0.75637579= 00 
0.762222A1E OC 
0 . 7677A597F 00 
0.773076*2E 00 
0.77835917E 00 
0. 73335 3A5E 00 
0 . 73 8 306C0E OC 
0.79322696E 00 
C. 79805553= OC 
0.8028A357E 00 
0 • 80 7596A5E 00 
0. 812315 3AE 00 
0.81695 6 16E 00 
C.9215A87CE 00 
0 . 3259201 OE 00 
0. 93020508E 00 
0.83422661= 00 
0. 83805382E 00 
0.8A1852A3E OC 
0.84538698 E 00 
0.843 39 174E 00 
0.R5237265F 00 
0.35577905E 00 
0.3591*685= 00 
0.8625C901E 00 
0. 95593793F 00 
0. 3691424 IF 00 
0.87241 943E OC 
0.37562978E 00 
0. 37379717c 00 
0.83175356E 00 
0.8 346502 8E 00 
0.S8729680E OC 
0.839851A5E 00 
0., 89232504= 00 
0.89^636=1= 00 
0.89693129= 00 
0.89920104= 00 
0.90141773= 00 
0.90361393E OC 
0. 9058G404E 00 
0.93797A84E 00 
0.9101 2955E OC 
0.91226876E 00 
0.91439247= 00 
0.91 648936E 00 
0.91354370= 00 
0.920526 39E 00 
0 • 92 24009 5E 00 
0. 92420399E 00 
0.92579663= 00 
0.92736032E 00 
0.92385 37 5E 00 
0.93026102= 00 
0.93165315c 00 
0.933026C8E 00 
0.534 3695 6E 00 
0. 93570352E 00 
0.93702793E 00 
0.938 3434 IE 00 
0.9396493AE 00 
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C . 952867365 00 
C. 354667345 00 
0.856667316 00 
C.35886729E 00 
0.860867265 00 
C. 862867246 00 
0.364867215 00 
C.86686718F 00 
0.868867165 00 
0.87C867135 00 
0. 8728671 IE 00 
C. 874867085 00 
0.876867065 00 
C. 878867Q3E 00 
0 . 88C86700E 00 
0.882866985 00 
0.884866555 00 
0 . 886866535 00 
G.38386S90E 00 
0.890866835 00 
C. 892866855 00 
0.894866825 00 
C. 8^6366806 00 
C. 308866775 00 
C. 900366755 CO 
0.902866725 00 
C.904S667CE 00 



0.940946345 00 
0.942234405 00 
0.943513515 GO 
0.944783095 00 
0.946039565 00 
0.947269205 00 
0.943463085 00 
0.949572326 OC 
0. °50663C9E 00 
0.951586375 00 
0 • 95 2 5048 75 OC 
0.953367355 OC 
0.954169635 00 
0.954967145 00 
0.955781545 00 
0.956515675 OC 
0.957274445 00 
0.958029036 00 
0.95977066c 00 
0.959525115 CO 
0.960265406 00 
0.961002116 00 
0.961734065 00 
0.962462435 OC 
0.963134245 00 
0.963906055 OC 
0.969618335 00 
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